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ABSTRACT 



Galactic globular clusters have been pivotal in our understanding of many astrophysical phenomena. Here we publish the extracted 
stellar parameters from a recent large spectroscopic survey of ten globular clusters. A brief review of the project is also presented. 
Stellar parameters have been extracted from individual stellar spectra using both a modified version of the Radial Velocity Experiment 
(RAVE) pipeline and a pipeline based on the parameter estimation method of RAVE. We publish here all parameters extracted from 
both pipelines. We calibrate the metallicity and convert this to [Fe/H] for each star and, furthermore, we compare the velocities and 
velocity dispersions of the Galactic stars in each field to the Besanfon Galaxy model. We find that the model does not correspond 
well with the data, indicating that the model is probably of little use for comparisons with pencil beam survey data such as this. 

Key words. (Galaxy:) globular clusters: individual: M4, M12, M22, M30, M53, M55, M68, NGC 288, NGC 6752, 47 Tuc - Galaxy: 
halo 



1. Introduction 

Globular clusters (GCs) are intriguing astronomical objects for 
many reasons. They are invariably found surrounding spiral 
and elliptical galaxies, have been used as tracers of galactic 
potentials (e.g. iKissler-Patig et all 119991 : iGebhard t & Thomasl 
1.2009) a nd can be employed to test gravitatio nal theo- 
ries (e.g. iLane et al.L l2009t ISoUima & NipotiL l2010h . Despite 
decades of detailed examinati on, it is still uncertain how 
GCs formed (Bekki & Chiba, l Om iLipsc v & Plavchan, 2004; 
iMashchenko & Sills, 2005; GriflFen et al., 2010; Lee et al., 2010) 
and their dark matter (DM) content is s till debated , altho ugh they 
are usually considered to be DM-poor (iLane et al.ll2009i and ref- 
erences therein). Recently, several studies revealed an interesting 
possible flattening in the velocity dispersions of several GCs at 
large radii, reminiscent of large elliptical ga laxies, a signature o f 
the kinematics of DM-dominated objects d Scarpa et al.L l2007h . 
These results called into question either the paucity of DM in 
GCs or our understanding of the nature of the gravitational in- 
teraction at low accelerations (below aq ~ lO"'" ms"^). 

We performed detailed kinematic studies on ten Galactic 
GCs, with the goal of calculating independent velocity disper- 
sion p rofiles to determine wheth er these results could be repli- 
cated dLane et al 1 12OO9I llQlOa'b, hereafter Papers I, II and III 
respectively). In short, we found that the flattening of the ve- 
locity dispersion profiles shown in previous studies were not re- 
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producible, despite two of our clusters being chosen to overlap 
with earlier investigations, namely M30 (Paper I) and NGC 288 
(Paper III). Our results indicated that neither DM, nor a mod- 
ification of gravitational theory, were required to reconcile the 
observed kinematics of our target GCs with current theory. In 
addition, we extended a recent metallicity calibration technique 
for open and globular clusters using the equivalent widths of the 
calcium t riplet lines and the horizontal branch m agnitude of the 
cluster (Cole et al.L l2004t |WaiTen & Colel l2009b . to the K band 
magnitude of the tip of the Red Giant Branch. This is several 
magnitudes brighter and can, therefore, be used for more distant 
clusters, and, more importantly, for GCs with hot, blue, hori- 
zontal branches whose stars do not exhibit strong calcium triplet 
lines (Section |2]l. Furthermore, a broad measure of the strength 
of the Galactic tidal field at various Galactocentric distances was 
made by comparing the velocity dispersions of the external and 
internal parts of the GCs. On the basis that the member stars at 
large radii are affected more by the external gravitational field, 
we showed that any flattening of the velocity dispersion profile 
could be attributed to this external field and did not require DM 
or modifications to gravity. We also argued that the lack of tidal 
heating signatures in our GCs at large Galactocentric radii is 
weakly suggestive of a spherical dark Halo (Paper III). In addi- 
tion, we discovered the exciting possibility of a two-component 
kinematic population within 47 Tucanae, which we interpreted 
as evidence for the cluster forming as two individual clumps in 
the protocluster cloud which coalesced at a later date (with an 
upper limit of ~ 7.3 + 1.5 Gyr ago; Lane et al., 2010c). We also 
found an as yet unexplained anomalously rapid cooling of the 
outer regions of GCs following tidal heating by the Galactic disc. 
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From these surveys we produced the largest sample to date 
of spectral data of both 47 Tucanae (47 Tuc) and M55 members, 
as well as large numbers of members of M4, Ml 2, M22, M30, 
M53, M68, NGC 288 and NGC 6752. Although the main goal of 
this manuscript is to make these data public for further research, 
from these surveys we also obtained spectra from Galactic field 
stars in the foreground of the GCs, therefore, we also have a 
large sample of Galactic spectra at various latitudes and longi- 
tudes. In the current paper we compare the velocity dispersions 
of the Galactic stars at each location to current dynamical models 
of the Milky Way. Projects such as HERME (Hi gh Res olution 
Multi-object Echelle Spectrograph; Barden et all |2008) on the 
Anglo-Australian Telescope show the importance of radial ve- 
locity/metallicity surveys over a large range of Galactic coor- 
dinates (particularly moving from the Disc to the Halo) for un- 
derstanding Galactic evolution. Here we present our pencil beam 
survey as an analysis of the Galactic velocity dispersion as a pre- 
cursor to surveys such as HERMES. Fuithermore, a comparison 
is made between metallicities calculated using analysis, as 
used by the Radial Velocity Experiment (RAVE; 'Zwitter et al.', 
|2008|). and the equivalent width method used by Lane et al. 
(I2010ah . Because our sample of stars from the 47 Tuc field is 
by far the largest, and also contains two GCs from the Small 
Magellanic Cloud (NGC 121 and Kron 3), we have used this 
field for most discussions and figures in the current paper. 

2. Derived Parameters 

All data was obtained on the AAOmega instrument on the 
Anglo-Australian Telescope (AAT), with the same gratings. We 
used the 1700D grating {R - 10000) on the red arm to obsei-ve 
the calcium triplet region (~ 8340A-8840A) and the 1500V 
grating (R - 3700) on the blue arm to include the swathe of 
iron and magnesium lines around ~ 5200A. A detailed descrip- 
tion of the observations, the membership selection process, and 
other details of the programme were given in Papers I, II and III. 

Our data pipeline (see Kiss et al., 2007), which calculates 
various stellar parameters, is based on the RAVE pipeline. All 
data from this programme was reduced using both the RAVE 
pipeline (which has been modified for use with AAOmega data) 
and our own for comparison. These two pipelines work in 
slightly different, albeit similar, manners. The Kiss et al. (2007) 
pipeline uses an iterative pr ocess to obtain best fits to synthetic 
spectra from the library by iMunari et al.l (l2005l) . which are de- 
graded to the resolution of AAOmega, and this model is cross- 
coiTelated with the observed spectra to calculate the stellar pa- 
rameters. The RAVE method obtains the best fit templates using 
penalised rather than cross correlation (a detailed description 
of the template matching, and subsequent parameter esti mation 
proce sses f or each are given by Kiss et al., 2007; Zwitter^tal], 
12008 ^). The ' Kiss et all (12007 ') version of the pipeline also differs 
slightly from the modified RAVE pipeline in that it trilinearly 
interpolates the spectra in the synthetic library (refining the grid 
in Teff, log g and [m/H]). This leads to resolution of 50 K in Teff , 
0.1 in log g and 0. 1 in [m/H], whereas the RAVE version reduces 
to a nonlinear interpolation on the six dimensi ons of the param- 
eter space. Fuithermore, The lKiss et al.l(l2007l) pipeline enforces 
one more iteration than the RAVE pipeline in the fitting pro- 
cess when calculating the radial velocities, which may provide a 
small improvment in the quality of the radial velocity estimate. 

In most cases the RAVE extracted parameters are very simi- 
lar to those from our own pipeline, however, there is a small but 

' ihttp : //www ■ aao . gov . au/AAO/HERMES/ 1 



significant offset between the uncalibrated metalUcity ([m/H]) 
extracted from the RAVE pipeline and those from ours (Figure 
[TJ. This is a known limitation of the RAVE pipeline which 
overestimates the metallicity for low metalliciti es due to the 
noise model assuming the same S/N for all pixels (IZwitter et al.L 
2008). Due to these subtle differences between the methods, and 
parameter estimates, we publish here all parameters from both 
pipelines. Published parameters are shown in Tables [T]|2] and [3j 
all data are available via the CDS. 

2.1. Parameter Uncertainties 

JHK magnitudes are taken directly from the 2 Micron All Sky 
Survey (2MA SS), and, therefore, h ave an uncertainty of ~ 0.03 
magnitudes ( Skrut skie et al. 112006'). V and I magnitudes are esti- 
mated from 2MASS JHK using unpublished transformations by 
G. Bakos (private communication). We have verified that these 
transformations have an uncertainty of V ~ 0.2 and I ~ Q . l mag - 
nitudes by cross correlation with data bv lWeldrake et al.l(l2004l) . 
who provide VI photometry of 43 067 stars within ~ 30' of 47 
Tuc to better than 0.03 magnitudes. 

Rotational velocity (vrot) estimates are theoreti cal because 
they c ome directly from the synthetic spectra by Munari et al.l 
(l2005h . however, the lack of any inclination information means 
Viot can be considered v-sin(/). We do not derive a formal uncer- 
tainy on the quoted values of Vrot, but because the two pipelines 
derive Vrot directly from the synthetic spectra, we can estimate 
an uncertainty by simple comparison between the two pipelines. 
The mean difference of v,ot between the two pipelines is ~ 
12kms"'. Taken at face value this can be considered the un- 
certainty on Vrot, however, due to the spectral resolution of the 
observations we advise that values of Viot ^ 30kms"' are much 
less reliable. 

Similarly, because Tgff and log g are taken directly from the 
synthetic spectra for both pipelines, and because we do not de- 
rive any formal uncertainties for these parameters, the mean dif- 
ferences between the two pipelines (276.5 K and 0.6, respec- 
tively, see Figure|3ll can be taken as the face value uncertainties. 

At the resolution of our observations, and due to the density 
of the template spectra, microturbulence values of < 2 km s"' are 
not resolved. Therefore we recommend that this be taken as the 
minimum uncertainty on the microturbulence values, although 
the true uncertainties are likely to be larger. 

The uncertainties in all metallicities ([m/H], [M/H] and 
[Fe/H]) are +0.1 dex (see Section l3T2l and Figure|5]l. 

Note that both [a/Fe] and microturbulence estimates are 
highly unrel i able a nd should be treated with caution (see 
Zwi tter et al.L 12008 , for additional details of parameter re- 
liability and uncertainties for parameter estimates from the 
RAVE pipeline). Fuithermore, the velocity uncertainties from 
the RAVE pipeline are the errors on the fits of the maxima of the 
correlation functions using a quadratic function. They are, there- 
fore, not physical uncertainties on the radial velocities and tend 
to overerestimate the true uncertainties by about 20%. The ra- 
dial velocity uncertainties from the Kiss pipeline are the formal 
errors from gaussian fits to the cross correlation function pro- 
files. Again, these are not physical uncertainties and may also 
slightly overestimate the true uncertainties, however, the diff'er- 
ence between the velocity estimates themselves from the diff'er- 
ent pipelines is small. Although the overall mean difference be- 
tween the two pipelines for the 47 Tuc field is ~ 4.2 km s"', this 
is mainly due to the large uncertainties in radial velocity esti- 
mates for the hot stars and reduces to ~ 0.3 km s"' when compar- 
ing only those stars selected as members (see Figure |2}; derived 
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Table 1. Parameters published in the current paper from the 
RAVE pipeline for all fields. The final column designates 
whether the star was classified as a member in Papers I, II and 
III. See text for explanations of the parameters and associated 
uncertainties. 



Cluster/Field 

RA (radians) 

dec (radians) 

estimated I magnitude 

field name 

fibre number 

V, (kms-') 

Vr uncertainty (kms"') 

^eff (K) 

[m/H] 
[ff/Fe] 

microturbulence (kms"') 
rotational velocity (v-sin[i]) 
cluster member (yes/no) 



Table 2. Parameters published in the current paper from the Kiss 
pipeUne (with fields centered on M22, M30, M53 and M68). 
The final column designates whether the star was classified as 
a member in Papers I, II and III. See text for explanations of the 
parameters and associated uncertainties. 

Cluster/Field 
Star ID 

V, (kms-i) 

Vr uncertainty (kms"') 

sum of the equivalent widths of the CaT lines (A) 
RA (degrees) 
dec (degrees) 
estimated V magnitude 

Teff (K) 

[m/H] 

rotational velocity (v-sin[i]) 
distance from cluster centre (") 
position angle 
cluster member (yes/no) 



velocities for stars with T^s > 9000 K are much less reliable as 
these have the CaT lines ove rtaken by strong PI 3 , PIS, and PI6 
hydrogen Paschen lines (e.g. lFremat et al.LII996l) . 

The RAVE collaboration represents the state of the art in the 
extraction of many parameters from stellar spectra. The small 
deviations between the parameters extracted with the modified 
RAVE pipeline and those extracted using our own pipeline show 
that our software is also representative at this level (Figures |2] 
andEJ. 



Table 3. Parameters published in the cuiTent paper from the Kiss 
pipehne (with fields centered on 47 Tuc, M12, M4, M55, NGC 
288, NGC 6752). The final column designates whether the star 
was classified as a member in Papers I, II and III. See text for 
explanations of the parameters and associated uncertainties. 

Cluster/Field 
Star ID 

V, (kms-') 

Vj uncertainty (km s"') 

Sum of the equivalent widths of the CaT lines (A) 

RA (degrees) 

dec (degrees) 

estimated I magnitude 

estimated V magnitude 

Teff (K) 

log,? 
[m/H] 

rotational velocity (v-sin[i]) 
J mag 
H mag 
K mag 

distance from cluster centre (") 

position angle 

cluster member (yes/no) 
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Fig. 1. Comparison of [m/H] from the RAVE pipeline versus that 
from our own for all stars in the 47 Tuc field. Notice the small 
offset between the parameter values. 



3. Results 

3.1. Kinematics of Galactic Field Stars 

Despite the main focus of this paper being the publication of the 
data, we include some analysis of the field stars here. To do this, 
it was first necessary to exclude all stars considered members of 
the GCs (i.e. extract the field stars from the complete dataset). 
Our strict membership selection method in Papers I, II and III 
evidently did not extract all cluster members because it is ob- 
vious some have been left behind (Figure |4]i. It was, therefore. 



necessary to impose an additional cut based on the distance from 
the centre of the cluster to ensure we removed all cluster mem- 
bers before the model comparisons. 

Comparing our results with those of the best available dy- 
namic Galaxy model (the Besangon model; iRobin et"an . I2003L 
see also http : //model . obs-besancon . £r/) gives a measure 
of the accuracy of the model, which is known to have limitations 
(e.g. the Disc component truncates at lOkpc, also see Figure 8 
by Conn et al., 2008). Table|4]shows the velocity dispersion and 
mean velocities of all non-members from our ten fields and the 
Besangon model. The model fields have the same field centres 
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Fig. 2. Comparison of radial velocity estimates from the RAVE 
pipeline versus that from our own for all the stars in the 47 Tuc 
field. Comparing only those stars determined to be members, the 
mean difference reduces to ~ 0.3 km s ' (Paper II). This is due 
to hot stars having strong Paschen lines (see text). 



as our data fields and were chosen to contain a similar number 
of stars, with the same temperature ranges as our data. The ve- 
locity dispersions for the data were calculated using the Markov 
Chain Monte Carlo method described in Paper III. The lack of 
agreement between the model and our data is apparent in Table 
m however, the Besangon model is regu larly employed for com 



pariso n with observational surveys (e.g. t onn et 



jioy ea tor com- 
dl l2005ll2007l 



120081) . The deviation of the model from the data in the current 
paper highlights the model's limitations for comparisons with 
pencil beam surveys. This should be taken into consideration 
when using this model for analysis. The Besangon model is, 
however, very eff' ective when used for comparisons with large 
area surveys (e.g. lRevle et al.Ll20Toh . 

A possible source of the model's discrepancies with small 
field surveys is that, while the model does include the warp and 
flare of the Disc, it does n ot include any known ti dal streams, of 
which there are many (e.g.'Belok urov etaP . '2007V Furthermore, 
note an apparent kinematically distinct population of stars with 
100 < V, < 200km s-i in the 47 Tuc field (Figure H. The 
Small MageUanic Cloud is in this field (Paper II) which also 
has a reces s ional velocity of 100 < Vr $ 200kms~' (e.g . 
Storm et al.|. l2004t lEvans & Howarthi 120081; iDe Propris et all 



,2010). The model cannot be expected to repUcate the mean ve- 
locity, or dispersion, of a field such as this. 

3.2. Metallicity Calibration Based on [m/H], Tgg and log g 

The RAVE project calibrated [M/H] and [Fe/H] b ased on [m/H1 
[a/Fe] and \ogg (see Equations 19, 20 and 21 by IZ witter et al.L 
l2008l who found that, using this method, the results for both 
[M/H] and [Fe/H] were within ~ 0.2 dex of reference metallic- 
ities). We have used this same method to calculate [M/H] and 
[Fe/H] for all stars in the current survey. Since our pipeline 
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Fig. 3. Comparisons between the estimated T^s and log g from 
each of the pipelines. Note that there is no obvious systematic 
offset between the two and that the average diff'erence is small. 



does not calculate [a/Fe] we cannot directly compare [Fe/H] 
calibrations based on our pipeline to those using the modified 
RAVE pipeline. However, [Fe/H] was calculated for the ten 
members of NGC 121 found in the 47 Tuc field in Paper II. 
These ten stars can also be seen in Figure |5] as the overden- 

sity at [Fe/H] 1.6. In Paper II the metallicity for these stars 

were calculated using an equivalent width calibration method 
which found [Fe/H]= -1.50 ± 0.10 for NGC 121, completely 
consistent with the overdensity in Figure |5] Furthermore, Figure 
|5]also shows the derived [Fe/H], using the RAVE methodology, 
for all stars from the M68 field. The overdensity centered on 
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Fig. 4. Stars observed within in the 2dF field centered on 47 Tuc. 
Red points indicate those stars selected as members in Paper II 
and black points indicate those selected as non-members. Note 
the overdensity of black points within the velocity range of the 
cluster which are likely to be members. Circled black points in- 
dicate stars selected as non-members for the current paper 

Table 4. Comaprison between the velocity dispersion (cr) and 
systemic (mean) velocity of non-cluster member stars in each 
field and the Besan§on model fields. All values are given in 
kms"'. 



Cluster/field 


0^ data 




Vsys{data) 


^sys(model) 


M4 


61.7 


78.2 


-15.1 


-36.2 


M12 


53.2 


68.1 


-20.4 


0.9 


M22 


60.9 


81.3 


-2.8 


32.6 


M30 


49.0 


52.6 


-10.9 


-9.4 


M53 


53.0 


40.0 


0.1 


-2.6 


M55 


59.4 


74.6 


-8.8 


7.6 


M68 


35.9 


51.4 


7.2 


12.8 


NGC 288 


40.8 


38.1 


7.8 


6.1 


NGC 6752 


50.4 


62.6 


-5.5 


-15.5 


47 Tuc 


58.8 


49.8 


39.1 


16.0 



[Fe/H] 2 are the stars determined to be members of M68 in 

Paper II, and are well separated in metallicity from the Galactic 
field stars. Using the equivalent width method, in Paper II we 
calculated [Fe/H]= -2.06 + 0.15 for M68, again c onsistent with 
that de rived in this paper with the RAVE method. ' Zwitter et al.l 
(I2008h stated that it was not possible (or at least very diflicult) 
to provide a physical explanation for the calibration relation be- 
tween metallicities derived from equivalent widths or photome- 
try methods to those obtained hy analysis. While this is still 
true, the current paper provides further evidence that the calibra- 
tion relations bv .Zwitter et al., (.2008.) do, in fact, hold. 

4. Conclusions 

We present for publication the stellar parameters of 29 35 1 stars, 
observed with AAOmega in fields centered on ten Galactic glob- 
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Fig. 5. The [Fe/H] and [M/H] calibrations for all stars in the 47 
Tuc and M68 fields. These match well with the values of [Fe/H] 
calculated in Papers II & III via the equivalent widths of the 
Call triplet lines and K magnitude of the Tip of the Red Giant 
Branch. The red circles represent [M/H] as the ordinate, and the 
black points represent [Fe/H] as the ordinate. A representative 
error bar is shown in the lower right of each panel. 



ular clusters. Parameters were extracted using two pipelines, 
namely a version of the RAVE pipeline adapted to work with 
AAOmega data and a pipeline based on the RAVE methodology, 
and we publish here all parameters from both pipelines. 

In addition, we find that the velocities, and velocity disper- 
sions, of our Galactic field stars (those not belonging to the clus- 
ters) in each field do not agree well with those of the Besangon 
Galaxy model. This descrepancy may be due to the model lack- 
ing information on the many tidal features present in the Galactic 
halo, as well as nearby objects like the Magellanic Clouds. It is 
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apparent that there is a population which is kinematically dis- 
tinct from the Galaxy in the 47 Tuc field (with 100 < V, < 
200km s '). Since the SMC is in the background of this field, 
which has a radial velocity of 100 < Vi < 200km s ', it is clear 
that we are seeing SMC stars in this field. We, therefore, suggest 
care is taken when using the Besan9on model for comparison 
with pencil beam surveys. 

Furthermore, we have calculated calibrated metallicities 
([M/H] an d rFe/H]) for e ach st ar based on the RAVE method 
outlined bv lZwitter et al.l (l2008l) and these correspond well with 
[Fe/H] calculated via the equivalent widths of the calcium triplet 
lines originally published in Papers I, II and III. 
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